Multicolor vortex solitons in two-dimensional photonic lattices 
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Abstract 

We report on the existence and stability of multicolor lattice vortex solitons constituted by cou- 
pled fundamental frequency and second-harmonic waves in optical lattices in quadratic nonlinear 
media. It is shown that the solitons are stable almost in the entire domain of their existence, and 
that the instability domain decreases with the increase of the lattice depth. We also show the 
generation of the solitons, and the feasibility of the concept of lattice soliton algebra. 

PACS numbers: 42.65.Tg, 42.65.Wi, 42.79.Gn 
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Localized structures, especially solitons, play a crucial role in many branches of the non- 
linear science. Over the past several decades, the existence and unique properties of spatial, 
temporal and spatiotemporal optical solitons in homogeneous cubic and quadratic nonlin- 
ear media have been studied both theoretically and experimentally (for detailed reviews, 
sec Solitons arise as a result of the balance between linear diffracting and/or 

dispersing properties of the medium and a nonlinear mechanism responsible for focusing 
/defocusing. One important subject of study is the generation of nonlinear modes with a 
nontrivial phase, such as vortices. In optics, vortices are associated with screw phase dis- 
locations nested in light beams |4]. Here we are interested in vortices with a bright shape, 
i.e., dislocations nested in finite-size beams 0. In homogeneous pure cubic and quadratic 
uniform media such ring-shaped vortex solitons suffer azimuthal instabilities which have 

n 

been observed experimentally in different settings p. They can be made stable in media 
with competing nonlinearities and in media with refractive index modulations that we 
address here. 

Propagation of optical radiation in media with transverse refractive index modulation 
differs considerably from the propagation in uniform media. Localized structures in such 
periodic media, termed discrete or lattice solitons, do exist and exhibit a rich variety of 
topologies. Since the theoretical prediction of discrete optical solitons in 1988 they keep 
attracting a growing interest, in part because of their potential for all-optical switching and 
routing phenomena l9j, llO, llll • The intermediate regime between continuous and discrete 
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solitons [12, 13] constituted by continuous nonlinear media with an imprinted transverse 
modulation of the refractive index, has been shown recently to offer a variety of new oppor- 
tunities. The concept behind this regime might be termed tunable discreteness, the strength 
of modulation being the parameter that tunes the system properties from continuous to 
discrete. In this context, wave dynamics is governed by the interplay between optical tun- 
nelling to adjacent sites and nonlinearity. This kind of lattice solitons have been observed 



recently in two-dimensional (2D) photorefractive optical lattices |l4i Il5l llfil ll7L Il8||. 

Besides the study of the existence and stability properties of fundamental (ground state) 
modes in the nonlinear media with periodic potentials an intriguing question is whether the 
action of the confining potential can permit the formation of vorticity carrying localized 
structures. Theoretical works showed that such complex localized structures, i.e., lattice 
vortex solitons exist when an optical lattice acts on a Kerr or photorefractive nonlinear 
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crystal 19Ll20l l2l|. Recently, theoretical expectations were indeed confirmed experimentally 
by two independent groups (3, [2^. During the last years various families of solitons in 
arrays of weakly cou pled waveguides made with quadratic nonlinear media have been also 



investigated 



24. 
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and observed for the first time recently [30(. One 



dimensional (ID) multicolor solitons in lattices with tunable strength have been also studied 
recently, and their potential applications for packing and steering single solitons have been 



investigated 
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Two-dimensional geometries might support even robust soliton ensembles 
with phase dislocations, a problem not analyzed so far, and which is the aim of this work. 

We thus report on the existence and stability of such multicolor lattice vortex solitons, 
which comprise four main humps arranged in a square configuration. It is shown that the 
lattice vortex solitons are stable almost in the entire domain of their existence except for a 
narrow region near the cutoff, and that the instability domain decreases with the increase 
of the lattice depth. We also investigate the possibility of their dynamical generation 3^1 
from Gaussian-type input beams with nested vorticities. 

We study the system of coupled nonlinear equations that describe the interaction between 
the fundamental frequency (FF) and second-harmonic (SH) waves under conditions for type 
I second-harmonic generation in bulk materials in the absence of Poynting vector walk-off: 



.dqi d x (d 2 q 1 d 2 qA , 

= — I + I - qiQ2 exp {-zp£) - pR (77, C) qi 



Mi 



2 \dr] 2 d( 2 
2 



+ 



- ql exp (iPO - 2pR (r), () q 2 , 



(1) 



<9£ 2 \dr] 2 d( 2 , 

where q\ = (2ki/k 2 ) 1 ^ 2 [2iru} 2 x^r 2 /c 2 ]Ai and q 2 = [2iruj 2 x^rQ/c 2 ]A2 represent the normal- 
ized complex amplitudes of the FF and SH fields, k\ = k(u)o), k 2 = k(2uJo) ~ 2k\, tq is the 
transverse scale of the input beams, r\ = x/r , ( = y/r Q , £ = z/(k\rl), (3 = {2ki — k2)k\r 2 is 
the phase mismatch, d\ = —1, d 2 = —ki/k 2 ~ —1/2, and p = 2injj 2 5x^ v> r 2 ./ 'c 2 is the lattice 
depth. The function R(rj,Q = cos(2irr]/T) cos(27r£/T) describes the transverse refractive 
index profile, where T is the modulation period. The system (1) admits several conserved 
quantities, including the energy flow 



and the Hamiltonian 



H 
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where V = e v (d/drj) + e^(d/d(), and e^, are unity vectors along 77 and £ axes. 

We searched for the stationary solutions in the form q± = (ui(r), () + ivi(r], ()) exp(i&i£) 
and g 2 = (1*2(77, C) + i^iVi 0) ex P(^20) where «i,2(?7, () and t>i,2(?7,C) are re& l functions, 
and 6i j2 are real propagation constants that verify 6 2 = ft + 2&i. Substitution of the above 
expressions into Eq. (1) yields the following system of equations for the soliton profiles n 12 
and t>i )2 

di / d 2 U! d 2 uA 
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We solved the system of coupled equations (4) numerically by using a standard re- 
laxation method. The lattice vortex soliton families are one-parameter families defined 
by the propagation constant h\ for any given period of the modulation T, lattice depth 
p and phase mismatch (3. Since one can use scaling transformations (71,2(77, Cj £, PiP) ~ > 
X 2 Qi,2(x r lj xC, X 2 £, X 2 Pi X 2 p) to obtain various families of solitons from a given family, we 
have selected the transverse scale r such that the modulation period is given by T = 7r/2 
and then we have varied b\, (3, and p. 

The simplest vortex soliton with unit topological charge in two-dimensional periodic lat- 
tice is shown in Fig. 1. It comprises four main humps arranged in a square configuration 
with a stair-like phase structure that is topologically equivalent to the phase of a conven- 
tional vortex in uniform medium (see Fig. l(b,d)). The positions of the soliton intensity 
maxima almost coincide with the positions of the local maxima of the lattice. Note that 
the singularity of these vortex solitons is centered between four lattice sites, that is they 
belong to the class of the, so called, off-site vortex soliton. In the model we investigated here 



there exist also a family of on-site vortex so 
singularity is centered on a lattice site 



itons (not shown here). In that case the phase 
We will restrict ourselves here to the 



x solitons (: 
0, 22, Q 



case of the off-site vortex solitons. It is interesting to note that these stationary structures 
somehow resembles the 4-soliton molecules carrying orbital angular momentum that were 
investigated in a variety of nonlinear media in both two-dimensional and three-dimensional 
34 , [35 1 . The typical stair-like phase distribution in the case of the above 



geometries 



33, 



mentioned soliton molecules is clearly seen in panels (b) and (d) of Fig. 1 for the lattice 
vortex solitons. We want to mention that the ID quadratic waveguides were also shown 
to support various families of multipeaked solitons, which display combinations of in-phase 
and out-of-phase odd solitons, the latter ones with n phase jumps between neighbor solitons 

0. 

It should be noted that at low powers (small b{) the lattice vortex solitons are quite 
wide and spread out over many lattice sites (Fig. 1(a)), while at high powers the energy is 
mainly localized within the corresponding four peaks (Fig. 1(c)). We did not find four-hump 
structures with higher topological charges (2 or more), and all other higher-order stationary 
structures we have found (for example, lattice vortex solitons with eight humps) were found 
to be unstable on propagation. Thus, in this work we will restrict ourselves to the study of 
the properties of simplest four-hump lattice vortex solitons. 

In order to characterize the families of lattice vortex solitons we have calculated the 
energy flows associated with these stationary solutions as well as their existence domains for 
given lattice depth p and phase mismatch (3. As a general rule, the energy flow of the four- 
hump lattice vortex solitons is a non-monotonic function of the propagation constant (even 
if this cannot be seen directly from Fig. 2(a) without zooming). Note that in Fig. 2(a,d) 
on the abscise we have plotted the difference (61 — b co ) between the propagation constant b\ 
and cutoff value b co . The cutoff value depends on both the phase mismatch (3 and lattice 
depth p. For example at p = 8 cutoff is given by b co = 6.065 for (3 = —6, while b co = 1.23 for 
(3 = 6. The cutoff b co is a non-monotonic function of the lattice depth p [Fig. 2(b)]. It tends 
to infinity at p — > and p — > +00. One can see from Fig. 2(c) that in the presence of the 
lattice the dependence b co {(3) differs from that for quadratic solitons in continuous media: 
bco(P) = max{— (3/2, 0}. Thus at (3 — > —00 cut-off is approximately given by ((3 — /?o)/2, 
while at /3 — > +00 one has b co = b Q , where (3q is the mismatch shift due to the lattice. Both 
60 and (3q growth with p. This property holds also for the one- dimensional lattice solitons 
in quadratic nonlinear media |3jJ. 

The periodic refractive index modulation affects also the energy sharing between FF and 
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SH waves. For example, at a given phase mismatch (3, the fraction of the total energy flow 
carried by the SH wave increases with increase of the lattice depth. Moreover, near the 
cut-off, the SH wave spreads over more lattice sites than the FF beam. As in the case of a 
uniform media, in the lattice with fixed depth p the part of energy flow carried by the SH 
wave decreases with increase of phase mismatch j3. 

To investigate the stability of the lattice vortex solitons, we have performed extensive 
numerical simulations of the evolution dictated by Eq. (1) with the input conditions q± (£ = 
0) = ( Wi (t7, C)+ivi(v, C))(l+Pifa 0) and g 2 (£ = 0) = (u 2 (r], ()+iv 2 (r], C))(l+p 2 (?7, 0), where 
ui t 2 and f i ; 2 are the exact solutions of Eq. (4) and pi j2 are random functions with Gaussian 
distribution and variance cr^ oise = 0.01. We have propagated the perturbed four-hump lattice 
vortex solitons over thousands of units for various values of the physical parameters involved 
(P, P, U). 

Our simulations show that there exists a narrow instability band near the propagation 
constant cutoff b co for vortex solitons, but above certain critical value of propagation constant 
they appear to become free of instability. We have found that the width of instability domain 
of lattice vortex solitons decreases with increase of the depth of the lattice p [Fig. 2(d)]. For 
example, as depicted in Fig. 2(d), the width of instability domain on propagation constant 
for p = 12 is approximately given by 0.24, while for p = 8 it is 0.35. 

A few representative decay scenarios for the unstable four-hump lattice vortex solitons 
with unit topological charge are shown in Fig. 3. In the row (a) of Fig. 3 we show the typical 
decay of the unstable vortex soliton in the vicinity of the cutoff on propagation constant 
b co . The initial energy of the localized structure is spreading out during evolution across the 
whole lattice and the vortex soliton disappears. Notice that this type of instability develops 
exponentially. In the rest part of the instability domain located closer to critical value of 
the propagation constant, we encountered oscillatory- type instability. Upon development of 
this instability vortex soliton transforms into a fundamental (ground state) lattice soliton 
that is the most robust and energetically stable state of the system (see row (b) of Fig. 3), 
through increasing oscillations of four intensity maxima of the vortex. 

One of the important results of this study is that the lattice vortex soliton becomes com- 
pletely stable when its propagation constant exceeds a critical value b cr , i.e., almost in the 
entire existence domain (see Fig. 2(d)). In row (c) of Fig. 3 we have plotted, for the sake of 
illustration, the initial and the final (after 500 propagation units) intensity distributions of a 
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stable lattice vortex soliton. Comparing to the soliton molecules investigated in bulk nonlin- 
ear media, which were shown to be metastable physical objects under suitable conditions, we 
conclude that, as expected on physical grounds, the effect of the two-dimensional lattice is 
to arrest the rotation of the soliton molecule and thus to assure the complete stabilization of 
the soliton complex. Since lattice causes strong azimuthal modulation of the vortex soliton, 
lattice removal results in complete soliton decay into four filaments, as shown in Fig. 4. 
Escape angles of filaments decrease with increase of input energy flow of vortex soliton. 

To understand lattice vortex solitons generation from a radially symmetric input beam 
carrying a screw phase dislocation nested in the center and to show that different sets 
of output solitons can be obtained with different combinations of topological charges and 
shapes of the input beams we performed a comprehensive set of simulations of Eq. (1) with 
the input conditions corresponding to Gaussian beams with a phase dislocation nested in 
the center: 

Qi(C = 0, r, = Ar' mi exp^mxy?) exp(— r 2 /u> 2 ) 

g 2 (£ = 0, r, <p) = Br^ exp(im 2 ip) exp(-r 2 /w 2 2 ), (5) 

where r = [rf + C 2 ) 1 ^ 2 is the radius, (p is the azimuthal angle, A and B are amplitudes of 
FF and SH waves , wi and w 2 are beam widths. Below we set the width wi — w 2 — 1 and 
suppose that topological charge of FF wave is given by m\ — 1. 

First, we consider non-seeded vortex soliton generation at B = 0. At low input powers 
both input FF wave and generated SH wave exhibit complete diffraction, and input beam 
energy is redistributed between many lattice sites. With the increase of input energy flow 
(i.e. by increasing A in Eqs. (5)) the generation of lattice vortex soliton with unit topological 
charge becomes possible as shown in Fig. 5. Lattice soliton generation is accompanied with 
energy radiation (Fig. 5(c) and (d)) but the ratio between radiative losses and the output 
soliton energy flow decreases with increase of input energy flow. 

In the case of seeded SH generation (B ^ and 5 <C A), the output field distribution 
can be controlled by the input topological charge of SH wave. For mi = 1 vortex soliton 
generation is possible only for the vorticity-matched case when m 2 = 2, while all other 
values of m 2 correspond to formation of trivial-phase soliton distributions, whose structure 
is dictated by lattice symmetry and energy exchange between FF and SH waves at the 
initial stage of propagation. Some representative output distributions are shown in Fig. 6. 
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These plot s show that the concept of " soliton algebra" previously explored in homogeneous 
media 36[ , does also apply in the presence of lattices, offering interesting opportunities for 
controlling the soliton dynamics. 

In summary, we have shown that periodic lattices imprinted in quadratic nonlinear media 
can support four-hump vortex solitons with unit topological charge that are stable provided 
that their propagation constant is above a certain critical value. Below this critical value 
we have identified two types of instabilities: (i) an exponential-type of instability leading to 
the final decay and spread out of the solitons across the lattice and (ii) an oscillatory-type 
instability leading to the transformation of the lattice vortex soliton into a fundamental 
soliton without internal vorticity. We investigated the generation of the multicolor lattice 
vortex soliton from Gaussian beams with nested phase dislocations. The possibility to 
generate different output lattice soliton patterns, with and without vorticity, by varying the 
topological charges and amplitudes of the input beams in seeded excitation configurations, 
has been discussed. The generation of a 2D periodic potential in quadratic nonlinear media 
is a challenging issue, even though fabrication of ID lattices has been already achieved using 
techniques which might be extended to 2D geometries. Also, the results presented here 
might be relevant to suitable atomic-molecular Bose-Einstein condensates held in optical 
lattices. 
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grant BFM2002-2861. 
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Figure Captions 

Fig. 1. (a) Profile and (b) phase of vortex solitons supported by the harmonic lattice at 
bi = 1.07. (c) Profile and (d) phase of vortex soliton at &i = 2. Only the FF wave is shown. 
Lattice depth p — 4, phase mismatch (3 — 0. 

Fig. 2. (a) Vortex soliton energy flow versus propagation constant for different values of 
phase mismatch at p — 8. (b) Propagation constant cutoff versus lattice depth at (3 — 0. (c) 
Cutoff versus phase mismatch at p — 8. (d) Stability and instability domains for different 
lattice depths at (3 — 0. Circles show critical value of propagation for stabilization. 

Fig. 3. Propagation of vortex solitons with b\ = 3.1 (a), 3.4 (b), and 5 (c) in the presence 
of input noise with variance o 2 noise = 0.01. FF wave profile is shown at different propagation 
distances. Lattice depth p — 8, phase mismatch (3 = 0. 

Fig. 4. Snap-shot images showing decay of the stable vortex solitons caused by removal 
of the lattice. Only SH wave profile is shown. Images are taken after each 2.5 propagation 
units. Lattice depth p = 8, phase mismatch (3 = 0. 

Fig. 5. Generation of the vortex solitons with only FF input, (a) Field and (b) phase 
distributions of the input FF beam with topological charge m\ = \. (c) FF beam and (d) 
SH beam at £ = 15. Lattice depth p = 8, phase mismatch (3 = 0. 

Fig. 6. Soliton algebra. The output soliton distribution depends on the topological 
charges mi of FF wave and m 2 of SH wave, respectively. In all cases, m x = 1. In (a)-(d), 
the amplitude of FF wave A = 20 and the amplitude of SH wave B = 2. In (e) and (f) 
A = 20 and B = 0.5. Plots (a)-(f) correspond the topological charges m 2 = 1,3,4,6,7,8 
respectively and show the output SH field distribution at £ = 100. Lattice depth p = 8, 
phase mismatch (3 = 0. 
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